"Stochastic Modeling of Coercivity " - A Measure of Non-equilibrium State. 
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A typical coercivity versus particle size curve for magnetic nanoparticles has been explained 
by using the Gilbert equation followed by the corresponding Fokker Plank equation. Kramer's 
treatment has been employed to explain the increase in coercivity in the single domain region. 
The single to multi-domain transformation has been assumed to explain the decrease in coercive 
field beyond a certain particle size. The justification for using Langevin theory of paramagnetism 
(including anisotropy energy) to fit the M vs H curve is discussed. The super-symmetric Hamiltonian 
approach is used to find out the relaxation time for the spins (making an angle greater than 90'^ 
with applied field) at domain wall. The main advantage of our technique is that we can easily take 
into account the time of measurement as we usually do in realistic measurement. 

PACS numbers: 75.60.Ej, 75. 75. -fa, 75.60.Lr. 



I. INTRODUCTION 



Thermal excitation and relaxation phenomena play a 
very crucial role in the case of nanoparticles. At finite 
temperature, single domain nanoparticles often exhibit 
superparamagnetic behavior, i.e., the relaxation time of 
the particles is much smaller than the characteristic time 
scale of the measuring instrument. Superparamagnetic 
behavior has currently been studied by a number of 
experimental techniques such as ac and dc susceptibility 
measurements 0, 0, H 0: neutron diffraction Q etc. 
These effects have considerable technological interest 
because of their relevance to the stability of information 
stored in the form of magnetized particles. 
Coercivity is an important quantity which plays a 
crucial role as far as the stabilization of a magnetic 
system is concerned. Understanding the nature of the 
coercive field with the variation of particle size is one 
of the central issues, which has been addressed in the 
present paper, by using stochastic theory 0. The 
behavior of the coercivity as a function of particle size 
is a old oneQ- The physical phenomena creating the 
maximum in coercivity of nanoparticles as a function 
of size is a well known problem in magnetism. This 
problem has been qualitatively understood. Various 
theoretical models have been published on the particle 
size dependence of coercivity |^ 1^ • Thermal switching 
in single-domain particles (where coercivity increases 
with size of particles) was considered by a lot of authors 
din in 111 El- These models failed to explain the de- 
crease in He with the increase in V at large particle size. 
Nucleation of domain walls was investigated by Braun 
|l5l Il6l |. The crossover from single to multi-domain 
switching was investigated numerically by Hinzke et al 
The effect of measurement time (i.e. the time lag 
between the measurement and the application of field) 
was not included in their approaches. In this paper 
we have emphasized to give a mathematical basis to 
explain this well known phenomena including the effect 
of measurement time. We quantitatively explain the 
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FIG. 1: Coercive force as a function of particle size. ( W. H. 
Meiklejohn, Rev. Mod. Phys. 25, 302(1953); F. E. Luborsky, 
J. Appl. Phys. 32, 1718(1961).) 



peak in the coercivity vs. particle size curve. 
In this paper we have tried to explain the non-monotonic 
(first increase and then decrease in coercivity, FIG.l) 
behavior in coercivity against the particle size with the 
help of non-equilibrium statistical mechanics approach. 
The time of measurement is automatically included 
in this approach. Our description is based on the 
assumption that our system is consist of mono-dispersed 
particles with no interparticle interaction. Although the 
particle size distribution and interparticle interaction 
can produce many interesting effects 0, ^ 01 1 but it 
is beyond the scope of this paper. 

With the preceding background the paper is organized 
as follows. In Sec. II we have discussed the increasing 
part of the He by assuming the material consist of 
single domain particle with high anisotropy barrier 
limit. This section also contains the magnetization 
calculation taking into account the anisotropy energy. 
In Sec. Ill we have discussed the decreasing part of the 
coercivity by using super symmetric quantum mechanics 
(SUSYQM) approach. Finally, in Sec.V we present our 
main conclusions about the significance of the reported 
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FIG. 2: Direction of the applied magnetic field, magnetic 
easy axis and magnetization. 
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FIG. 3: Energy profile of an uniaxial single domain mag- 
netic particle (a) in absence and (b) in presence of external 
magnetic field. 



II. SINGLE DOMAIN REGIME 

The subject of how a bulk magnetic specimen acquires 
a single domain structure and exhibits magnetic viscos- 
ity due to Neel relaxation, when its size is reduced, is an 
old one ■ The critical volume of a single domain par- 
ticle with uniaxial high anisotropy energy as estimated 
by Kittel is i?c = (9cr„/47r/5^), where (t„ is the do- 
main wall energy per unit area and Is is the saturation 
magnetization of the particle. The numerical value of 
Rc comes out to be ^ 10 nm for K ^ 10^ erg/cc (K is 
the anisotropy energy per unit volume). Therefore the 
initial increase in cocrcivity is observed in single domain 
particles. This leads us to map our problem to that of 
an ensemble of uniaxial single domain particles in con- 
tact with a heat bath at temperature T. Let us assume 
that the magnetic moment vector of each particle (/2p) 
(FIG. 2) makes an angle 6m with the easy axis and this 
easy axis of the particle is at an angle 9k with respect 
to the apphed magnetic field(^app)- The total energy of 
the system is given by: 



Et = KVsir?{9m) - fipHappCos{ek - 9„ 



(1) 



The Gilbert equation governing the dynamics of the 
spin is 



dt 



dflp 



V- 



dM 
'dt 



(2) 



The energy profile of the above system (for Happ=Q) 
is a symmetric double well potential (with two minima 
at 6*1 and 92, generically referred to as 9i and a max- 
ima at 03(FIG. 3)). Here we shall apply high barrier 
and weak noise limit to simplify our mathematical calcu- 
lations. Hence Kramer's ansatz can be used to find the 
probability of getting a particle making an angle 9m with 
the easy axis and is given by p3 | 



P{9m,t) ^ P{9,,t)e- 



(3) 



Although we can solve the corresponding Fokker Plank 
equation with the help of SUSYQM. approach, we have 
used an approximate method to get a better feeling of 
the underlying physics of the problem. The number of 
particles oriented between 9i — 9m to 9i + 9m is given by 



ni = P{9i,t)e "bt 



where, 



e *'bt sin{9m)d9r, 



(4) 



(5) 



with i = 1 or 2. Since in our case the potential barrier 
{5E) is much higher than the thermal fluctuation, we 
have used the saddle point approximation in Eq. (5) which 
results in 



e ^^ksT 



(6) 



where Ci is the curvature at the i-th minima. Now the 
continuity equation for P{9m,t) is given by 



dP{9m,t) 
dt 



-V.J', 



where, 



J' = J - kVP{9m,t). 



(7) 



(8) 



The second term of the above equation represent a dif- 
fusion in 9m — P{9m) space, which takes care of thermal 

agitation and J~vP{9m,t), where v = We have 

also assumed the diffusion coefficient to be space inde- 
pendent. In the limit of high relaxation time J'(0„i,t) 
becomes independent of 9 and is given by 



sin{9m) 



1 9ET{9m)p^^^^^^^ , dP{9m,t) 



ksT d9r, 



d9r, 



(9) 
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FIG. 4: Computer simulated magnetization curve of two dif- 
ferent particles (a)7 nm (b) 8.5 nm, using Eq. (11). 



Using Kramer's ansatz in the above equation one can 
write 



dn2 _ 

where = ToI:iIi and I3 
tion of Eq.(lO) is 



n2 _ 

hi 



(10) 



-dOm- The solu- 



n2{t) = 



nT2 



-e '■"i'''^2)*(nT2 - n2(T2 -l-ri)) 



Tl + T2 



(11) 



where (rii — ri2) is proportional to the net magnetiza- 
tion along the direction of the applied magnetic field. 
Eq.(ll) needs one initial condition i.e. 712(0) to get the 
value of n2{t). For a single domain particle with large 
relaxation time, if one changes the magnetic field after 



a finite interval of time (t). then limgH- 



H~5H 



(0) ^ limsH- 



it) 



(0). This implies 



that for a particular value of H one should not expect 
to get the same value of magnetization during increas- 
ing and decreasing cycle of H. Since the relaxation time 
Ti increases with particle volume, — also in- 

creases with particle volume giving rise to higher coerciv- 
ity. Hence coercivity is a consequence of the quasi-static 
non-equilibrium measurement. Therefore, Langevin the- 
ory of paramagnetism is not applicable in these cases. 
We have used Eq.(ll) to numerically generate M vs H 
curve as shown in FIG. 4, for particle size 7 nm and 8.5 
nm. We have used t ~ 50scc and K=10^erg/cc in the nu- 
merical calculation, which is realistic for measurements 
of coercivity by vibrating sample magnetometer. 
In case of a superparamagnetic system, a common prac- 
tice is to fit the magnetization curve by using Langevin 
theory of paramagnetism |25l | . But it is not quite obvious, 
since Langevin theory of paramagnetism does not include 



the magnetocrystaline anisotropy energy. Therefore an 
attempt has been made here to calculate the magneti- 
zation of a superparamagnetic sample at equilibriums. 
Magnetization 



< M >= fip < cos9 >, 
where, < cos9 J^^ P{9, Happ) cosO sinO d9, 

P(9,Happ) - "2 I 

Et = —Hjip COS0 + k sin^9, 
and partition function 



(12) 



Z = 



,-^^(i^./.[VA:,(- + i)] 



(13) 



where k — KV (V is the volume of the particle), a ~ 
Happ£p_ i,^ _ _k_^ Here Erfi is the error function 

k-Bi 



and fci 

26j. The magnetization of the system 



< M >= ^J.pS[kl,a] - — , 



with 



S[ki,a] = ^J,p- 



X 



(14) 



(15) 



where X ^/k;^{Erfi[Vk,{§: + 1)] - Erfi[Vki{§: - 1)]). 
Now < M > ^ fipL{a) at the limit ki — > 0, where 
L(a) is the Langevin function. Therefore in general it 
is better to use Eq.(14) rather than Langevin theory of 
paramagnetism. But for very small particle size one can 
use Langevin theory of paramagnetism directly. 
One can calculate the magnetization for the superpara- 
magnetic state using the rate equation at equilibrium 
for the particles whose barrier height is at least 5 times 
larger than the thermal energy. Let us consider the 
situation where we have applied a dc magnetic field 
parallel {9k=0) to the easy axis. If 71,2 and rii be the 
number of particles having spins parallel {9k = 0) 
and antiparallel {9k = 180°) to the applied field, the 
magnetization along the direction of magnetic field is 
Ms{n2 — ni), Ms being the saturation magnetization. 
From the master equation at equilibriums we have 



< M{n2-ni) >= Ms{ 



T2-T1 



where Ei = E, 
above equation gives us 
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< M, >= RL 



(16) 
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FIG. 5: Spin distribution near 180'' domain wall (a) at real 
space (b) at probability space. 



This indicates that the anisotropy energy does not affect 
the magnetization of superparamagnetic particle (with 
large anisotropy energy) in equilibrium. 



III. MULTI-DOMAIN REGIME - 
SUPERSYMETRIC QUANTUM MECHANICS 
(SUSYQM) APPROACH 

Let us now concentrate on the second region where 
the coercivity field decreases with the increase in parti- 
cle size. It is clear from the above discussion that this 
can not happen if the particles still comprise single do- 
mains. Since the coercivity of the single domain particle 
increases monotonically with the increase in particle vol- 
ume, hence a single to multi domain transformation takes 
place at the maximum of coercivity. 
To explain the decreasing part of the coercivity field let 
us consider an arrangement of spin in a linear chain as 
shown in FIG. 5. In the following discussion one should 
keep in mind that we arc not interested here in the origin 
of the domain wall, but we assume the existence of do- 
main and the Hamiltonian contains the relevant terms. 
The Gilbert equation corresponding to i-th spin is 



~dt 



dt 



(18) 



Let us consider cylindrical polar coordinates, where z = 
represents 0th site, z = a represents 1st site and so on. 
For a particular site the spin can move over the surface of 
the cylinder along a semicircular curve {Oe[Q, tt]). Writing 
everything in terms of cylindrical polar coordinate and 
proceeding as above we have 



di 



with h' 



70 ntJ- 



1 dVjd) 



subjected to the 



do 
(19) 

condition 



r(e,t) 

dz 



0, 



1+70 J?/^' 

which physically means that there is no 
spin hopping between two sites. This implies that the 
spin will start relaxing along the surface of the cylinder 
without changing its position along Z-axis. To solve the 



above equation let us perform the following transforma- 
tion P{d,t) = y/l%(W)2P{9,t). Where P^gie) = Aqb'^ 
and e ~ ksT. Using the transformation in Eq.(19) we 



have ^ = + 



dt 

function 



v"{e) 

2 



^j— ) ip. Introducing a new 



such that ip{6, t) 



17 



, we obtain 



(20) 



with ^ ^ defining two operator A 

= —-§g+^, such that eA^Acf) = Xcj) with ground state 
eigenvalue equal to zero (since A(j)Q = to get equilibrium 
distribution). It can be shown that if (pi be the first 
excited eigenstate of A^A then it is the ground state of 
AA^ with ground state eigenvalue Ai (23|. Now one can 
apply variational method to get Ai 



d_ 

de 



^ and 



Ai = 



/0i(6')e^At0^(g)rf6i 



(21) 



Now we are in a position to get the solution of the 
Fokker Planck equation. Suppose we have applied a 
static magnetic field along the crystallographic direction 
of the sample. The potential energy of the spins making 
an angle less than 90*^ can be approximated by harmonic 
oscillator like potential ■^Kd^. Which gives us the 
relaxation time of the order of hi{e — K). The spins 
at the end of the domain wall have a double well like 
potential. For such double well potential we have 
Ai - ;ii(e-(^°-^(''i» + e-(^<'-^(«2)), where Vo is the 
barrier height and 9i and 02 are the position of the 
two minima. Now it is clear that the relaxation time 
depends on the damping parameter as well as the barrier 
height, which in turn depends on the value of anisotropy 
constant and the angle between successive spins. The 
anisotropy constant is higher (one order of magnitude) 
for smaller particle. So for a smaller particle it takes 
more time to relax back to its initial configuration, 
giving rise to a higher coercive field. The above model 
indicates that our system under consideration, consist 
of a linear chain of ferromagnetic particles having two 
domains, with their easy axes parallel to each other 
and with applied magnetic field also. The above model 
also does not contain the domain of closure. Still the 
model can be regarded as the starting point to explain 
qualitatively the hysteresis of a multi domain system. 



IV. CONCLUSION 

The particle size dependence of the coercive field is 
explained from the view point of non-equilibrium statis- 
tical mechanics. The zero hysteresis loss indicates that 
the system reaches its equilibrium state before the next 
data is taken. It is shown that the increase in coercivity 
is the effect of increase in relaxation time. The numer- 
ical calculation of hysteresis curve shows that the loss 
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as well as the coercive field increases with particle size. 
During the numerical calculation assumption has been 
made that ksT < SE. The Langevin theory of paramag- 
netism for superparamagnetic particle has been reestab- 
lished, after taking into account the magnetocrystaline 
anisotropy effect. Assuming the single to multi-domain 
transformation we have shown that the relaxation time 
of the samples decreases with increases in particle size 
due to decrease in surface pressure and anisotropy con- 
stant which gives a decrease in coercivity. The important 
point to note that in case of any experimental study of a 
single domain particle, one should perform this coercive 
field (He) vs. particle size curve and figure out the peak 
in the curve. One should perform all the measurements 
below this peak particle size to analyze the behavior of 



single domain magnetic nanoparticles. 
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